A nonlinear quantum piston for the controlled generation of vortex rings 
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We propose a simple way to manage interactions in Bose-Einstein condensates to generate vortex 
rings in a highly controllable way. The vortex rings are generated under the action of a quantum 
analogue of a classical piston pushing the condensed atoms through a small aperture. 
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Introduction.- One of the most remarkable achieve- 
ments in quantum physics in the last decade was that of 
Bose-Einstein condensation in ultracold alkaline atomic 
gases. These physical systems have a high potential for 
supporting quantum nonlinear coherent excitations and 
in fact, many types of nonlinear waves have been experi- 
mentally observed [l" 12| or theoretically proposed to ex- 
ist (see e.g. the reviews [l3|)) in ultracold quantum de- 
generate gases, mainly with Bose-Einstein condensates 
(BECs). In this way BECs are, apart from their funda- 
mental interest, one of the best physical systems for the 
manifestation of nonlinear phenomena and certainly the 
best one in the quantum world. 

But in the jungle of nonlinear excitations the king 
is probably the vortex ring. Vortex rings are essen- 
tially three-dimensional topological nontrivial structures 
appearing in either classical or quantum 15[ flu- 
ids. They are able to propagate in cylindrically trapped 
BECs as stable objects (16] . similarly to what happens 
in classical superfluids [17|. This is an essential differ- 
ence with most other solitonic structures that become 
unstable when passing to fully three-dimensional scenar- 
ios, e.g. one dimensional bright solitons that are unsta- 
ble to blow-up f?] or dark solitons, that are unstable to 
the snake instability [1]. This leaves the vortex ring as 
the only dynamically nontrivial nonlinear excitation ob- 
served in BECs in three-dimensional scenarios. 

Vortex rings were first observed in BECs as the out- 
come of the decay of dark solitons [8|[and as a result of 
the decay of quantum shock waves [9|. More recently 
they have been observed to appear durin g co mplex oscil- 
lations in soliton-vortex ring structures |10| and during 
the merging of BEC condensate fragments However, 
all of those scenarios involve complicated nonlinear phe- 
nomena and in general a simple mechanism allowing for 
the controlled generation of a prescribed finite number of 
vortex rings is still missing. The purpose of this paper is 
to propose a mechanism allowing for the generation of a 
few vortex rings in a highly controllable way. 

Physical idea.- The process of vortex ring generation 



in classical fluids has received a substantial treatment in 
the literature. One of the most standard ways to obtain 
vortex rings in classical fluids involves moving a piston 
through a tube, resulting in a vortex ring being gener- 
ated at the tube exit. A standard generation geometry 
consists of the tube exit mounted flush with a wall with 
the piston stroke ending at the tube exit [19]. 

Is it possible to export this idea to the quantum world 
in order to obtain vortex rings in BECs? One might try 
to imagine a complicated setup to create a tube-like ge- 
ometry using ordinary potentials and having a moving 
part in order to generate the flow, but this would lead to 
a complicated trap geometry and would be difficult to ob- 
tain experimentally. In this paper we will use something 
conceptually much simpler using the possibilities open by 
space-dependent Feschbach resonance management in a 
Bose-Einstein condensate. Since the first achievements in 
scattering length control in BECs [l^, the technique of 
Feschbach resonance management has been improved and 
used in many different applications. Presently, the level 
of control of the scattering length allows for its very pre- 
cise tuning pl| and nothing prevents an extended control 
of the interactions leading to a space dependent scatter- 
ing length. A large number of theoretical papers have 
recently studied nonlinear phenomena in systems with 
managed interactions IJ, |22h24 | . 

The physical idea is very simple. Starting from an 
equilibrium BEC in a trap with a given value of the scat- 
tering length ai we propose modifying interactions in half 
of the space (say z < 0) to a larger value 02 > ai. This 
change would affect the initial configuration by inducing 
the transverse expansion of the atom cloud for 2; < 
while at the same time generating a fiow towards the re- 
gions with smaller interaction values located at z > 0. 
Those combined effects generate around z = a flow to- 
wards z > analogous to the piston-driven flow through 
an aperture used to generate vortex rings in classical flu- 
ids. It is interesting to note that all those effects are 
achieved simultaneously with a single action on the in- 
teractions without restorting to complicated external po- 
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FIG. 1: Atom flow when increasing interactions for z < 
((72) in an equilibrium condensate with gi — 10^, — 
Xy = l,\z = 0.05 according to Eq. ([1]). Shown are curves 
for (/2/ffi =1.6 (sohd brown), 1.8 (sohd black), 2 (solid 
red), 2.1 (dashed green), 2.2 (dashed blue). Shown are (a) 
the mass imbalance Mit{t) — J^^g\il){x,y,z,t)\^dxdydz — 
J^^^^\ip(x,y, z,t)\'^dxdydz as a function of time and (b) the 
flow Jz^Q(t) = pudS where 5 is a surface orthogonal to 
the unit vector u — (0, 0, u) at 2 = 0. 



tentials. In what follows we will show that this idea leads 
to a highly controllable method for the generation of a 
small number of vortex rings. 

Mathematical model.- We will consider a BEC in the 
mean field limit ruled by the mean-field equations 

2i^ = -Aij + V,^t{x)ij + g{zM^i;, (1) 

where the spatial coordinates and time t are measured 
in units of oq = h/2muj and l/2a;, respectively, while 
the energies and frequencies are measured in units o{2huj 
and 2a; respectively, uj being the transverse frequency of 
the potential in physical units and Vext = (A^a;^ + A^y^ + 
A^z^)/2. We assume that the total density is normalized 
to 1, so / iT/ipdT^ = 1. Also, g{z) = Anas{z)/ao is propor- 
tional to the local value of the s-wave scattering length 
0^(2:) and will be taken, starting from t — 0^, as 



5W = 




92 > 91- 



(2) 



and phase (f) [ip = p 



Analytical estimates. - 
Eq. ^ for density p 

consists of the continuity equation pt -I- V • (pu) 
and the integrated Bernoulli-type equation 2(/)( -I- 



The hydrodynamical form of 

) 




Kxt + g{z)p = V^^/^, where velocity is given by 
u = V0. These equations can be obtained from vari- 
ational principle <5 / / LdVdt — for the Lagrangian 
L = p{d^,cj)f/2 + (9..p)V8p + pV,xt + 9{z)p'/^ lead- 
ing to the momentum equation dt{pui) + dx^^ij — 0, 
where the momentum flux density tensor S^j is given by 
= (p + pd.j:iVcxt)Sij + pUiUj - pdl.^^ \tlp/A: and the 
pressure p = g{z)p'^/A. The local speed of sound is given 
by = dp/2dp = g{z)p/2. If t < 0, then g{z) = gi, 
and the system is in its ground state characterised by 
the chemical potential p, such that idtip = P'<P- For 
t < in the Thomas-Fermi approximation we have pQ 
(15giAiA2A3/1287r)2/5 and density po = (2/io - Kxt)/^!- 
As t — > 00 the system will reach a new ground state with 
Poo = (15gi52AiA2A3/647r(5i -|- 92))^^^, so that the den- 
sity will become 



Vcxt)/gi, 

Vcxt)/92, 



z>0 
z <0, 



(3) 



with 9i/{gi + g2) particles on the left and 92/(91 + 92) 
particles on the right of the plane z = 0. At time t = 0+ 
when the interaction strengths are turned to step func- 
tion given by Eq. ^ , we have equal number of particles 
on both sides of z = plane, so (52 — 5i)/2(5i + 92) par- 
ticles have to flow from the left half-plane to the right 
half plane before the equilibrium is reached. 

We consider the equations of conservation of mass and 
momentum in the neighbourhood of 2; = 0. We can 
rewrite the mass conservation in integral form as 



(4) 



Taking the limits — ^ gives piUi = P2U2, where index 
1 (2) characterizes velocities and densities at z — > 0+ 
(z — )> 0~). Similarly from the conservation of momentum 
equation and dropping the anisotropic quantum stress 
tensor we have piu\-\-gip\/ 4: = p2U2+g2P2/^- Combining 
these two equations gives 



P2 giPi 



92pI 



Pi 9iPi 



92pI 



Api pi - P2 



4/92 Pi - P2 



(5) 



Let's assume that ^ — 92/91 — 1 is small, so that 
(5=1 — p\/ P2 is also small. The requirement Ui — ^ 
as ^ — >■ leads to (5 = ^/2 and U2 = \J g\p\(,/2. From 
this it follows that U2 (note that U2 > ui) may exceed 
the speed of sound c if ^ > 1. So we expect the vor- 
tices to nucleate once 52/51 > 2. Our numerical simula- 
tions below confirm this prediction. The rate of change 
of mass across the plane z = can than be estimated 
as dM/dt = it2P2 ds where S2 is the cross-section at 
z — > 0~. Approximating pi by poo gives 



dM 



V(g2 -gi)gi 

dt 4 " (.g2 + .gi) 



(6) 
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FIG. 2: (Color online) Pseudocolor plots of the atom density 
\tp{r,z,t)\'^ for: (a) t = 0.15, (b) t = 1.5, (c) t = 2.25 and (d) 
t — 4. for gi/gi ~ 2.1. The spatial region shown corresponds 
to z £ [— 30, 30],r e [0,7]. Blue corresponds to low atom 
densities and red to high ones. 



Results.- To test our ideas we have run many numer- 
ical simulations of our model equations ^ for different 
parameter combinations. We have used a solver based on 
a fourth order finite difference scheme in space together 
with a fourth order Runge-Kutta discretization in time. 

In what follows we describe the typical outcomes for a 
specific example corresponding to a large repulsive BEC 
with gi — 10^, with = Xy = l,Xz = 0.05 (i.e. soft 
longitudinal trapping) using a computational window of 
240 space units in z and 10 units along the transverse di- 
rections x,y. Typical space steps are about 0.25 and the 
time steps At ~ 1.5 x 10"''. Our initial configuration is a 
ground state BEC corresponding to g{z) = gi. On that 
initial configuration we suddenly increase interactions for 
z < and observe the subsequent condensate dynamics. 

The results of a series of simulations to be described in 
detail later are summarized in Fig. [TJ Once the nonequi- 
librium situation is generated there is a flow of atoms 
from z < to z > depending on the values of the ratio 
92/91- For ratios 92/91 S [1,2] we find that our analyt- 
ical predictions for mass rate of change and velocities of 
the flow are quite accurate below criticality (eq. from 
Eq. ^ we get dM/dt ~ 0.01 for 52 ~ 1-6 as compared 
with Fig. [2]). However, when a predicted critical value 
92/91 — 2 is surpassed there are differences in the curves, 
such as a saturation in the mass transfer and a drastic 
reduction in the flow for certain times corresponding to 
the emission of vortex rings for those times as we will 
see later. Snapshots of the time evolution of the atom 




FIG. 3: (Color online) Same as in Fig. [3] but for (72/51 = 2.2 
and times (a) t = 0.15, (b) t = 16.5, (c) t = 21, (d) t = 33. 

density are shown in Figs. [Hill 

First, in Fig. [5] computed for 92/91 — 2.1 just around 
the critical value we observe that the initial spatial inho- 
mogeneities in the interactions leads to a steady atom 
flow directed towards the region with smaller interac- 
tions. For small times the sudden increase of interac- 
tions results in a strong flow of atoms coming from the 
region with z < towards the region with z > (right) 
resulting in the formation of a shock wave [23|. Shortly 
after that, a small vortex ring appears in the low density 
regions around z = but it counterflows towards z < 
and disappears in that low density unstable region. Af- 
ter that, the situation remains stationary with time with 
some additional depletion of the condensate density for 
z < due to the remnant flow directed towards z > 0. 

In Fig. [3] we present some results corresponding to 
92/91 — 2.2 just above the critical value for vortex ring 
formation. After the initial shock wave [Fig. EJa)] dis- 
appears from the central region, a vortex ring enters the 
condensate around z = coming from the low density 
regions corresponding to large r values [Fig. ^h)] and 
moves slowly through the condensate remaining stable for 
long times [Fig[3l[c,d)]. Another vortex ring with a large 
radius seems to be present in the lower density regions 
where it would be experimentally difficult to detect. 

Increasing the interactions even further to 92 /gi = 2.4 
leads to a richer dynamics as summarized in Fig. |4l The 
short-time dynamics is analogous to the previous cases 
[Fig. [3l[a)] but then a complex transient appears where 
several vortex rings enter the condensate and also rar- 
efaction pulses are clearly identified [see Fig. IH^b)]. After 
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FIG. 4: (Color online) Pseudocolor plots of the atom density 
\'ip{r,z,t)f for: (a) t = 0.15, (b) t = 7.5, (c) t = 20.25, (d) 
t = 26.25 and (e) t = 44.25 for 02/91 = 2.4. The spatial region 
shown corresponds to z £ [—30, 30], r £ [0, 7] in subplots (a-d) 
and to z £ [—10,50] in subplot (e). Blue corresponds to low 
atom densities and red to high ones. 

that, some of those vortices counter-flow and disappear 
and a very regular picture arises with several vortex rings 
moving to the right. Fig. Sl^c) shows two vortex rings 
slowly moving through the condensate for t — 20.25 and 
a third one being generated around z = 0. Fig. SJ^d) 
shows a later stage of the evolution where three long-lived 
vortex rings travel smoothly through the condensate al- 
though their relative positions change due to small differ- 
ences in their speeds (notice the small differences in their 
radii) and their interaction with sound waves originating 
after the reflection of the shock wave in the condensate 
boundary [see Fig. IH^e)]. 

Conclusions.- We have proposed a simple method to 
create a quantum piston able to generate vortex rings in 
Bose-Einstein condensates in a highly controllable way. 
Our proposal is accessible to present experimental tech- 
niques and improves essentially currently used methods 
to produce vortex rings that are based on the instability 
of nonlinear unstable structures and lead to a complex 
tangle of vortex lines and rings or to combined vortex 
rings and soliton complexes. 
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